Dynamics of a single vortex line in a Bose-Einstein condensate 
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We study experimentally the line of a single quantized vortex in a rotating prolate Bose-Einstein 
condensate confined by a harmonic potential. In agreement with predictions, we find that the vortex 
line is in most cases curved at the ends. We monitor the vortex line leaving the condensate. Its 
length is measured as a function of time and temperature. For a low temperature, the survival time 
can be as large as 10 seconds. The length of the line and its deviation from the center of the trap 
are related to the angular momentum per particle along the condensate axis. 
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The macroscopic motions of quantum and classical flu- 
ids are dramatically different. The description of a quan- 
tum fluid by a macroscopic wave function ip = yfp e %e 
imposes strong constraints upon its velocity field. At a 
position of non-vanishing density p, the velocity is re- 
lated to the phase by v = hVQ/m, where m is the 
mass of a particle of the fluid. Hence V x v = 0. A rota- 
tional motion of the fluid can be obtained only through 
the nucleation of vortex lines, along which the density is 
zero and around which the circulation of the velocity is 
quantized in units of h/m ||. 

Quantized vortices play an essential role in the dy- 
namics of all quantum macroscopic objects. Examples 
are flux lines in superconductors H, and vortex lines in 
superfluid liquid helium and gaseous Bose-Einstein 
condensates (BEC) §, |, [J |[ §. Among the remaining 
problems, is the shape of a vortex/flux line and the study 
of its time evolution. The observation of inclined flux 
lines in superconductors was only possible due to recent 
advances in electron microscopy [^0| . In gaseous BEC, 
one has comparably easy access to the vortex line because 
the density of the atom cloud is low. A few disordered 
vortex lines have been observed by taking tomographic 
images perpendicular to the long axis of a cigar shaped 
condensate J5J. An array of many vortex lines in a pan- 
cake shaped condensate has been observed by transverse 
imaging of the whole atom cloud pj| . 

In this Letter we report the full length observation of a 
single vortex line in a cigar shaped condensate. We find 
that as a result of spontaneous symmetry breaking the 
line is generally bent. Our experimental results confirm 
recent predictions, in which the shape of the vortex line 
minimizing the energy of the gas was derived for a given 
rotation frequency [H| Q [l^, jl6| . We also study the 
time evolution of the shape of the line. As the angular 
momentum of the gas slowly decays, the bending of the 
vortex line and its deviation from to the center of the trap 
increase. We investigate the influence of temperature 
on these dynamics. From our results one can hope to 
draw some indications for the shape and dynamics of the 
vortex/flux line in systems where a direct observation is 
not yet possible. 



Our 87 Rb condensate is formed by radio-frequency 
(rf) evaporation of 10 9 atoms in a Ioffe-Pritchard mag- 
netic trap. The atoms are spin-polarized in the F = 
ttif — 2 state. The magnetic trap has a longitudinal 
frequency uj z /2tt = 11.8 Hz and a transverse frequency 
(u x + u y )/4:7r = uj±/2tt = 97.3 Hz (the x axis is verti- 
cal). Because gravity slightly displaces the center of the 
trap with respect to the magnetic field minimum, the po- 
tential in the xy plane is not perfectly isotropic and we 
measure a 1% relative difference between uo x and uu y . 

The initial temperature of the cloud, which is pre- 
cooled using optical molasses, is 100 jaK. The conden- 
sation threshold is reached at T c ~ 300 nK, with 7V C ~ 
2 10 6 atoms. We cool to typically T ~ 90 nK which 
corresponds to a condensate with No ~ 5 10 5 atoms and 
a chemical potential /i ~ 70 nK. It is obtained using 

= ^o + 10 kHz as the final rf, where is the frequency 
at which the trap is emptied. During the rest of the ex- 
perimental cycle, we maintain the evaporation rf at an 
adjustable level, typically v = vq + 12 kHz. This allows 
us to control the temperature while observing the vortex 
line. 

Once the condensate is formed, we use an off-resonant 
laser beam to impose on the trapping potential an el- 
liptic anisotropy in the xy plane [||. The wavelength of 
the beam is 852 nm, its power 0.1 mW and its waist 
20 /im. Acousto-optic modulators deflect the position of 
the beam in the xy plane, thereby rotating the potential 
anisotropy at a frequency of Q/2tt = 70 Hz. We apply 
this "laser stirrer" for 300 ms, during which 7 or more 
vortices are nucleated |l7|| . After the stirring phase, the 
condensate evolves freely in the magnetic trap for an ad- 
justable time r. 

The preparation of a single vortex line takes advantage 
of the slight static anisotropy of our magnetic trap, so 
that the angular momentum is not exactly a constant of 
motion. In a time r ~ 1-2 s, we observe a transition from 
a multi- vortex condensate to a condensate with a single 
vortex. This relatively long time levels the fluctuations 
that may occur during the nucleation process. Thereby 
we are able to reproduce a single vortex condensate on 
every experimental cycle. This vortex line can then be 
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FIG. 1: (a) Schematic of the imaging system. The cigar 
shaped condensate is imaged after 25 ms of time of flight 
leading to the inversion of the ellipticity in the xz plane. Two 
beams image the atom cloud simultaneously along the longi- 
tudinal (z) and transverse (y) directions of the initial cigar, 
(b-d) Simultaneous longitudinal and transverse images of con- 
densates after r = 4 s (b), r = 7.5 s (c) and r = 5 s (d). 



studied for a time r < 10 s. 

The atom distribution at r is probed (destructively) 
by switching off the magnetic trap, letting the cloud ex- 
pand during £tof = 25 ms and performing absorption 
imaging. Two imaging beams aligned along the y and z 
directions probe the atom cloud simultaneously (Fig. [l]a). 
The beams are combined onto a camera with the same 
magnification. During the expansion, the transverse di- 



mensions x and y of the condensate are magnified by 
cjx^tof ~ 15, while the longitudinal dimension is nearly 
unchanged |l8| . It has been shown theoretically that the 
presence of a single vortex line does not alter this expan- 
sion and that the coordinates of the line are scaled by 
the same factors [ijfl . 

Figs, p] b-d show three condensates imaged after var- 
ious times. The left column shows the "longitudinal" 
view along z, representing the atom distribution in the 
xy plane. The right column depicts the "transverse" view 
taken along the y direction, representing the atom dis- 
tribution in the xz plane. The vertical (x) direction is 
identical for all images. The transverse images show the 
typical atom distribution where the ellipticity is inverted 
with respect to the in situ cigar form, caused by the trans- 
verse expansion during £tof- 

As in we use the longitudinal images in Fig. [l] to 
verify the presence of a single vortex. The transverse 
image in Fig. [l|b obtained for r = 4 s shows the vortex 
line as a line of lower atom density. Clearly this vortex 
line is not straight. It rather has the shape of a wide 
"U". The x position of the axial part of the vortex line 
(bottom of the U) is close to the center of the condensate. 
It coincides with the x position of the dip in density seen 
in the longitudinal image. In the laboratory frame, the 
curved vortex line is expected to rotate around the z axis 
with a frequency related to the angular momentum of the 
condensate. This is confirmed by the fact that we observe 
up- and downwards bending with equal probability. 

Fig. Qc shows images taken after r = 7.5 s, for which 
the angular momentum has decreased significantly com- 
pared to Fig. [I]b. In the longitudinal view, one sees a 
vortex off-center and in the transverse view a narrow U, 
where the bottom of the U no longer extends to the center 
of the condensate. 

Fig. |l|d shows a vortex line in the shape of an unfolded 
"N" , observed after r = 5 s. The width of the N, which 
is the projection of the vortex line onto the z axis, is 
comparable to the width of the U in Fig. |l|b. The fact 
that U as well as N shaped vortex lines are observed leads 
to the question whether the bending occurs in one plane 
or whether three dimensional deformations of the vortex 
line can also occur. We have indications that this may be 
the case: some transverse views (not shown here) reveal 
asymmetric U or N shaped vortex lines. 

In order to give a quantitative analysis, we measure 
by hand the distance along the z direction between the 
two points where the vortex line leaves the condensate 
(Fig. |^a). We do not distinguish U and N shaped vortices. 
Normalization by the length of the condensate along z 
leads to the quantity d z . In Fig. ^|b, we plot d z as a func- 
tion of r. Each point in the plot corresponds to a sin- 
gle image. The small spread of the points demonstrates 
the reproducibility of the vortex shape. The normalized 
length of the vortex line d z decreases quasi-linearly with 
time. 
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FIG. 2: (a) Schematic of the extraction of d z and d± from 
the vortex line, (b) Evolution of d z with r. Each point cor- 
responds to a single image. 
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FIG. 3: Variation of d z with the temperature T. 
from fits to data similar to those in Fig. gb. 
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Fig. |]b shows that after r = 7 s the vortex is still 
present. This differs from our earlier reports || simply 
because we have used two different methods of identifying 
the vortex. In |J only the longitudinal view was avail- 
able. Now we find that a vortex can still be identified in 
the transverse view, while in the longitudinal image the 
density variation due to the bent vortex line is within 
the background fluctuations. This long lifetime can be 
compared with the MIT result 0, where a vortex array 
with more than 100 vortices was produced at r = 0. The 
number of vortices was divided by 4 in ~ 5 s; however 
a single vortex could still be detected after r = 40 s. In 
both experiments, it is clear that the decay time of the 
last vortex is much longer than the decay time of the 
initial vortex lattice. 

We have repeated the experiment of Fig. || for different 
evaporation radiofrequencies v (implying different tem- 
peratures T) during the free evolution time r. We always 
observe a quasi-linear decrease of d z and fit its slope, 
d z . It is plotted in Fig. | as a function of T. Clearly 
the longest vortex lifetime is obtained for small tempera- 
tures, confirming the idea that the vortex line is dragged 
to the edge of the condensate by a thermally activated 
process Q . Note that at the lowest temperature we ob- 
serve almost exclusively U shaped vortices, while U and 
N are equally probable at the second lowest temperature. 




FIG. 4: Variation of d z and d± as functions of the average 
angular momentum per particle L z . All measurements were 
binned into intervals of L z = 0.1 H and averaged. The error 
bars give the statistical spread. 



In a last experiment we vary again r in order to prepare 
condensates with various lengths of the vortex line. We 
measure the angular momentum per particle L z by the 
same method as in |20| . The presence of a vortex leads 
to the lift of degeneracy between the two quadrupole 
surface modes carrying angular momentum ±2h 
00+ — ou- =2 (L z ) / (m (r_L 2 ) ) , where r± is the radius in the 
transverse plane. We excite a superposition of these two 
modes by a 1.9 ms flash of our laser stirrer and we probe 
the quadrupole oscillation at t osc = 2 ms and t osc = 9 ms 
after the flash. We repeat the experiment a third time 
without the laser flash and analyze the shape of the vor- 
tex line. As above we extract the normalized length d z of 
the vortex line along the z axis. We also measure the dis- 
placement d± of the axial part of the vortex line (bottom 
of the U) from the center of the condensate, and nor- 
malize it by the radius of the condensate in the xy plane 
(Fig. ||a). Since we have access only to the projection 
of the decentering on the xz plane, we actually measure 
d± | cos a\ , where a is the azimuthal angle of the axial part 
of the line. We account for this geometrical factor by di- 
viding the measured displacement by (| cosa|) = 2/tt. 

Fig. H shows d z and d± as functions of L z /H. For clar- 
ity we group all points into bins of L z =0.1 h and average 
over d z and d±. The error bars give the statistical vari- 
ation. The data corresponding to L z < 0.2 h are not 
reproducible enough and are omitted. The graph shows 
that d z and L z /h are approximately equal. A straight 
(d z ~ 1) and well centered (d_\_ <C 1) vortex line corre- 
sponds to an angular momentum of the order of h. When 
the angular momentum decreases, we measure a decen- 
tering d± < 0.15 as long as L z > 0.5 h. Below L z = 0.5 ft, 
d± rises to 0.3. 

We now compare our experimental results with recent 
predictions for the shape of a vortex line in an inhomoge- 
neous cigar shaped condensate |l|, |l4j [l|, |l6| . These 
theoretical studies consist in looking for the ground state 
of the condensate in a frame rotating at an angular fre- 
quency Q. The general conclusion is that above a critical 
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frequency Q c the ground state of the system has one or 
several vortices. The central vortex is generally bent if 
the trap aspect ratio oo±/co z is large compared to 1, which 
is the case in our experiment. A simple physical picture 
of the bending is given in ]l5[ . A cigar shaped conden- 
sate can be viewed as a series of 2D sheets of condensate 
at various z. For each sheet, one has a 2D vortex prob- 
lem leading to a critical frequency fic 2D ^ (z) above which 
a centered vortex is the stable solution. For a given ro- 
tation frequency it can happen that a centered vortex 
is the minimal energy configuration for the sheets close 
to z = (i.e. ^c 2D ^(0) < f)), while it is not the case 
for the sheets close to the edges of the condensate where 
the atom density is lower. In this case, the vortex line 
minimizing the total energy is well centered for \z\ < z c , 
and is strongly bent for \z\ > z c , where Q^ B \z c ) = 0. 
A precursor of this bending effect can also been found in 
|^2| , p3| in which a stability analysis of a straight vortex 
in an elongated condensate showed that some bending 
modes have negative energy, and are thus unstable. 

Our experimental procedure is somewhat different 
from the one considered in these theoretical studies. In 
our case no rotating anisotropy is imposed onto the con- 
densate during the relevant evolution. The stirring laser 
has been switched on for a short time only, at the begin- 
ning of the procedure, in order to set a non-zero angular 
momentum in the system. We observe the evolution of 
the condensate in our static trap, as the angular momen- 
tum of the gas slowly decays. However, the shape of 
the vortex line that we observe at short time in Fig. 0b 
is remarkably similar to those predicted and plotted in 
@ @ 111 111- As pointed out in [yj, this bending is 
a symmetry breaking effect which does not depend on 
the presence of a rotating anisotropy and which happens 
even in a completely symmetric setup [[24] . A decentered 
single vortex similar to the one shown in Fig. 0c was also 
found as the ground state of the rotating system in Q| 
for a given rotation range. By contrast, we did not find 
in the literature predictions for N-shaped vortices, such 
as the one shown in Fig. |l|d. This probably means that 
a N-shaped vortex is slightly more energetic than a U 
shaped vortex with the same angular momentum, and 
that it could not emerge from a procedure aiming to find 
the ground state of the system. 

Our results also provide information on the dynamics 
of the vortex line and the way it escapes from the conden- 
sate. A theoretical model has been proposed, in which 
the decay is due to the coupling with the non-rotating 
thermal component Jl9[ ] . In this model, the authors con- 
sidered a spatially homogenous condensate in a cylin- 
drical vessel. The vortex line, assumed to be straight, 
spirals out of the condensate. The bending of the line 
that we observe experimentally may change qualitatively 
the picture, since we find that the decay occurs first by 



an increased bending and, only in a second step, by a de- 
viation of the center of the line (bottom of the U) from 
the center of the condensate. 

In conclusion we have reported the observation of the 
full line of a single quantized vortex. We have related 
its shape (bending and deviation from center) to the an- 
gular momentum of the system. Our results should help 
modelling the dissipative evolution of a rotating Bose- 
Einstein condensate. 
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